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Four and Six Fermion Event Generators for e+e Collider Physics 

Elena Accomando 
Phys. Dept. Texas A & M, College Station 
TX 77843, USA 

The status of four and six fermion event generators for Standard Model processes 
at present and future e+e~ colliders is briefly reviewed. 

1 Introduction 

The ongoing run at LEP2 has led in the last rpjcent years to detailed analyses re- 
lated to both SM and beyond the SM physicsEl. The theoretical results achieved 
for LEP2 physics have been implemented in dedicated four fermion codes □. 
Mw measurements, studies of triple gauge couplings and Higgs searches in 
the low mass range require the analysis of four fermion final states, making 
therefore unavoidable the calculation of four fermion processes at LEP2 and 
future e+e" colliders. Since the advent of LEP2, several codes with different 
features have been implemented and carefully cross-checked. The available 
four fermion tree level programs, most part now interfaced to QCD Genera- 
tors (PYTHIA/JETSET and/or HERWIG), reviewed by Sjostrand, to relate 
parton level predictions with experiments, are listed in Tab.[l]. 

The work done for four fermion physics at LEP2 can be extended to the 
NLC. This machine however would not only improve the sensitivity to trilinear 
gauge couplings and gauge boson properties, but it would be also ideal for 
precision studies of quartic gaugg couplings, top properties and Higgs searches 
in the intermediate mass range El. All these topics involve processes with six 
fermions in the final state, but despite of this only a few six ferirdon codes have 
been at present implemented. To our knowledge, only GRACHD, SIXPHACTQ 
and WWGENPV/ALPHicI can simulate six fermion processes. 

In this talk we summarize some features of the above mentioned codes, 
even if this brief description can do no justice to the effort that has been done 
by the various groups. 

2 Four fermion codes 

On the side of four fermion codes, many improvements have been recentlv 
introduced into the existing programs (e.g. refer to contribiitions by JadachQ 
and Perret-Gallix) and two new codes, presented by Dennefl and Peskin, have 
been implemented. 
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Oode 


processes 


nif 


AG 


ALPHA 


all 


Y 


N 


CompHEP 


all 


Y 


Y 


ERATO 


CC11/CC20 


N 


Y 


EXCALIBUR 


all 


Y/N 


Y 


GENTLE/4fan 


CCll/NC 


Y/N 


Y 


grc4f 


all 


Y 


Y 


KORALW 


all 


Y 


Y 


LEPWW 


CC3/NC2 


N 


Y 


PANDORA 


CC3/NC2 




N 


PYTHIA6/JETSET 


New added 


Y/N 


N 


WPHACT 


all 


Y 


Y 


WTO 


all 


N 


N 


WWGENPV/HIGGSPV 


all 


Y/N 


N 


WWHP 


CC3 




Y 



Table 1: Available four fermions programs. The included processes are written using the 
notation of Ref. [1]. Y{N) indicates the code includes (does not) fermion masses (rrtj) and 
anomalous couplings (AC). Y/N is for approximate treatment of mf. 



Due to the increasing statistics at LEP2, the interplay between theory and 
experiments is becoming more stringent. The attention is therefore focusing 
on the role of the radiative effects and on the estimate of the theoretical uncer- 
tainty on the predictions. An example is the treatment of the photons in the 
initial and final state, a particularly delicate issue. Neglecting the px of the 
photon, as in the usual LL approximation of the ISR, can have sizeable effects, 
for example on the detection efficiency and on the differential distributions 
used for TGG studies and also for New Physics searches since the detection 
of new particles often rely upon missing px- Different prescriptions have been 
adopted. Some codes include models for the generation of IS multiphotons 
with finite Pt (KORALwIi) and LL FSR, while others make use of QED par- 
ton showers with pr (EXCALIBUR, grc4f, LEPWW and PYTHIA). A more 
accurate approach, not affected by possible gauge invariance problems, could 
consist in merging the explicit emission of a single hard photon with coUinear 
bremsstrahlung and radiative corrections. As a further step toward the full 
calculation of the radiative effects, complete matrix elements for the process 
e+£^ ^ 4/ + 7 have been implemented by ALPHA and recently by Denner et 
all. 

Matching theory and experiments means also being able to produce reliable 
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results in a finite time. This question has focused the attention on the inte- 
gration methods. Since a large number of diagrams with different topologies 
contribute to a given final state, the multichannel techniques (as in EXCAL- 
IBUR, KORALW and WWGENPV/ ALPHA) seem to be the most promising. 

At the higher NLC energies, not only the complexity of the processes will 
increase, but new difficulties could appear. In Fig.[l], it is shown the result of 
one of the latest tuned comparisons among four codes, concerning the gauge 
invariance issue for the single W production at small angle, reviewed by BoosU. 
Events with one or two electrons in the forward direction (coming from evW, 



500 GeV - No ISR 
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: cosa„>0.997 



: cose„>0.997 



2 3 4 5 6 



2 3 4 5 6 



Figure 1: Comparison for single W production at small angle at 500 GeV with and without 
ISR. On the X-axis: l=CompHEP, 2=grc4f, 3=KORALW, 4=WPHACT. Cuts: M{ud) > 5 
GeV and E^u, d) > 3 GeV. 



eeZ or 27 processes) will be relevant at the NLC as signal, to disantagle the 
WW J coupling, and as background to New Physics searches. These processes, 
which contain t-channel photon diagrams with mass singularity in the very 
forward region, require massive matrix elements and phase space to cure these 
apparent divergencies. Moreover, they need a proper implementation of the 
boson width in the propagators, in order to preserve the large cancellations 
among the t-channel diagrams dictated by the gauge invariance. The com- 
parison in Fig.[l] involves gauge preserving prescriptions, namely the overall 
prescription (CompHEP and WPHACT) and the L^^iy transform method (grc4f 
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and KORALW), and shows a rather satisfactory agreement. Recently, also a 
new scheme based on complex gauge boson masses and obeying the Ward 
identities has been introduced, as presented by DennerEl. 



3 Six fermion codes 

Six fermion final states receive contributions from a great amount of different 
Feynman diagrams. Some of them correspond to Higgs production and decay. 
Others may come from tt or WWZ, or from partial resonant processes (as for 
instance v^VeWW) as well as non resonant ones. All these diagrams, which 
in general are of the order of several hundreds, correspond to the same final 
state and interfere among themselves. Different resonant contributions might 
be considered as signal or background depending on the process under study 
and their interplay can be consistently analyzed only by complete calculations. 
Moreover, all these channels might constitute a main source of background to 
New Physics searches. 

Recently three groupiHl have started computing full tree level cross sec- 
tions for six fermion final states at the NI£!. TheseLXomputations have been 
applied to study the phenomenology of ttffl, WWz^ and higgffl production 
and have already shown the importance of finite width effects and irreducible 



background. In these analyses, the beamstrahlung e 
using the parametrization implemented in CIRCE 



sets have been simulated 
Recently, also a n£w 



method (PYBMS), based on the Pisin-Chcn formalism, has been proposedllil. 



4 Conclusions 

Much effort has already been devoted in order to have reliable tools for four 
fermion physics, but much work has still to be done on the side of electroweak 
and QCD radiative effects. Whereas the four fermion area is well covered, so 
far a few codes have been implemented for six fermion processes. Hopefully, 
the joint work of the various groups will lead to a better understanding of what 
is still missing in order to fully exploit the great and unique capability of the 
e+e^ colliders. 
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